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Holographic Fermi arcs and a d-wave gap 
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We study fermion correlators in a holographic superfluid with a d-wave (spin two) order parameter. 
We find that, with a suitable bulk Majorana coupling, the Fermi surface is anisotropically gapped. 
At low temperatures the gap shrinks to four nodal points. At high temperatures the Fermi surface 
is partially gapped generating four Fermi arcs. 
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Introduction and Summary. - Gauge gravity duaUty 
[IHS] provides a useful tool for understanding the strongly 
coupled planar limit of gauge theories. At low tempera- 
tures and non-zero density these theories exhibit proper- 
ties similar to the ones observed in condensed matter sys- 
tems. Using gauge gravity duality, various authors have 
constructed configurations which exhibit superconductiv- 
ity 111 H] , a p-wave order parameter [BHH] , and nodes in 
fermion correlators [9l[T0]. For reviews see e.g. [TltiTS] . 

In this work we pursue the analogy with condensed 
matter systems further and study fermion correlators in 
large A'' gauge theories with a d-wave (spin two) order pa- 
rameter. Our main finding is that in the phase where the 
spin two field condenses, and at non-zero temperature, 
the spectral function for the fermions develops asymmet- 
ric features including an angle dependent gap along the 
Fermi surface (FS). These features are exhibited in fig.jl] 
For particular values of the coupling between the spin 
two field and the fermions, this gap vanishes along a fi- 
nite interval reminiscent of the "Fermi arcs" observed in 
the pseudogap region of certain high superconductors 
[TH - flB] . We will discuss similarities and differences be- 
tween our results and the experimental observations in 
high Tc materials towards the end of this work. 

Setup. - Gauge gravity duality relates the partition 
function of a gravitational theory with an asymptotically 
anti-de-Sitter (AdS) solution to the partition function of 
a conformal field theory (CFT) in one spatial dimension 
less. We will refer to the gravitational theory as the bulk 
theory and to its CFT dual as the boundary theory. 



The bulk fields of a holographic superconductor with 
a d-wave order parameter in 3 -1- 1 dimensions are the 
gravitational field, a U{1) field and a massive charged 
(symmetric) spin two field ip^i, [17H19j . The unique mat- 
ter Lagrangian quadratic in i^^^, containing operators of 
dimension four or less, and describing the correct number 
of propagating degrees of freedom is [201 HI] 
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together with the constraint that the metric is Einstein. 
We introduced the notation (/?^ = D'^ip,^^ and (p = ip';^. 
The U{1) gauge field is denoted A and its field strength 
is F = dA. The Ricci scalar and Riemann tensor are R 
and R^^pa respectively. The fully covariant derivative 
is related to the metric covariant derivative through 

The Lagrangian ([T]) should be treated using the lan- 
guage of effective field theory. The equations of motion 
(EOMs) following from C exhibit non-causal behavior 
that can be corrected by adding higher order terms in 
F and its gradients. See [TH] for a detailed discussion. 

The Einstein condition on the metric restricts us to the 
probe limit, where the charge of the spin two field is large 
enough so that the matter fields do not backreact on the 
metric. Thus, the solution to the Einstein equations is 
the AdS4 Schwarzschild black hole whose line element is 
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FIG. 1: (Color online) A density plot of the fermion spectral 
function p defined in ( 14 1 evaluated at tj = for temperatures 
T = 0.49 7; (left) and T = 0.59 (right). Red and blue 
correspond to large and small values of the spectral function. 



ds-" 



(-/ dt^ + dx^ + dy^ + dz^/f) L'^/z^ (2) 



where f{z) = 1 — {4ttTz/3) and T is the Hawking tem- 
perature of the black hole and also the temperature of the 
field theory. To solve for the matter fields we consider an 
ansatz where only 



'Vvv = '^a(2) (3) 



are non-zero and (/j^^ is real. We note that with this 
ansatz, ip = ipp_ = F^^pLpP = 0. Under a spatial rota- 
tion of an angle in the (a;, j/)-plane, the spin two field 
transforms as 
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In [IH] it was noted that by redefining if^^ = z^ipi^^^ and 
choosing an angle 6 at which either (^a or (fxy vanishes, 
the EOMs and boundary conditions (BCs) for At and <p 
reduce to those of the AbeUan Higgs theory first discussed 
in [31 [S]. This identification guarantees that when the 
chemical potential /i of the boundary theory docs not 
vanish, there exists a critical temperature below which 
the spin two field condenses. 

We are after the spectral function p(aj, k) for ferniion 
operators O in the condensed phase of the boundary the- 
ory. A study of thermal fermion correlators in gauge 
gravity duality was initiated in and we will fol- 

low their approach. Consider the bulk fermion action 

+ r,*ipl^WT^D''^ - r|^T^'D''{ip^,^') , (5) 

where now — + \w^..,\a^^ ^ iqA^ with w the 
spin connection and \1/ a Dirac spinor. Gauge invariance 
of the coupling between the fermions and the spin two 
field dictates that the charge of ^' is q = q^/'i,. The 
bulk gamma matrices T'^ satisfy {r'^.T'^} = 25^''. The 
conjugate spinor '^'^ is defined by = CT-^*, where 
CT'^C^^ — — r^^. Vielbein indices are underlined. The 
phase of the coupling constant 77 can be absorbed in a 
redefinition of 5': we will use this freedom to set 77 > 0. 

If we restrict ourselves to couplings of mass dimension 
smaller than six then with the ansatz ([s]) and the choice 
q^ = 2q, the only other possible coupling between the 
spin two field and the fermions is (ci -I- 02^^) 4*. 

The latter provides an effective mass for the fermions 
and does not contribute to the asymmetric features of 
the spectral function which we would like to exhibit. 
Therefore we will not consider it further. Note that 
had one chosen q = qip, an interaction term of the form 
'pI.^'P^"* (c3 + caT^) \1/ would have been possible, lead- 
ing to an s-wave instead of a d-wave gap. 

The Dirac equation following from ([s]) is given by 

= (r'^D^ -m)^ + 2ir]ip^^T^'D''^'' + i-qLpf.V^" . (6) 

We find it convenient to introduce a rescaled spinor ■0 = 
{—g ■ g^^y^^'^. Since the Dirac equation couples ^ and 
\I''^, we consider a solution of the form 

iP = e-*"*+»^-^ C{z) + e*"*-^'=-^ C(z) , (7) 

where C and C can be decomposed into two two- 
component spinors: C = (Ci:C2)^- To write equation 
^ explicitly we must choose a particular representation 
for our spinors. Our bulk gamma matrices are given by 

pt ( —10-2 \ -pz f (J3 \ 

^ ~ V j<T2 y ~ V y ,n\ 

T^X_ f (71 \ T^y / — j(T2 ^ ' 

^ — \ ai) ^ \ia2 ) ■ 

In this representation C = T-. The boundary gamma 
matrices are defined by the action of F— on the positive 



eigenspace of F-: 7* = —ia2, 7^ = ci, 7^ = —0-3. Given 
([s]), the Dirac equation (|6| reduces to 

= £'(2) C2 + 2T]{g''''')hx[ipA<JlC2 + iVxy(72Cl] 

where we have set ky = 0. Arbitrary k can be recovered 
using The differential operators are defined as 

D{a) = V¥^<^3dz - m 

+ (-l)"(c^ + <zA,)v/^a2 + ifc:,Vff^ai . (10) 

The equations for can be obtained from (|9| with the 
substitution (w, kx) (— w, —k^) and C — > C- 

Once we find a solution to the Dirac equation in the 
bulk, we can use it to determine correlators of fermion op- 
erators O in the boundary theory. In particular, we will 
be interested in the Fourier transform of the retarded 
Green's function Gi?(i,f) = ie(t)({0(t, x), 0^(0)}), 
from which we can compute the spectral function. In 
what follows we summarize the procedure for obtaining 
the boundary theory Green's function. The reader inter- 
ested in a detailed exposition is referred to |22) . 

Consider the near boundary (small z) series expansion 
of the fermion fields, 

»°(ow)''"'+((.%1>"'"' <"> 

where the index a = 1, 2 is a spinor index. The compo- 
nent Sot acts as a source term for the fermion boundary 
operator Oa and Ra = (Oa)- A similar expression in- 
volving Ra and Sa follows from a series expansion of 
C- Since the equations of motion are linear, the two- 
component spinor R is linearly related to 5* and S |36| . 

Ra^MjSp+M^^S; . (12) 

The prescription for computing the retarded Green's 
fimction Gr{ui, k) is 

Gr = -iMj' (13) 

where the EOMs ^ should be solved with BCs which 
are infalling at the horizon. With a slight abuse of ter- 
minology, we define the spectral function p as 

p = TrImG'fl (14) 

where the trace is over the spin degrees of freedom. 

The appearance of a gap. - To compute p we need 
to solve the Dirac equation. However, even without an 
explicit solution to ^ we can infer the appearance of an 
angle dependent gap in p whose size vanishes along four 
nodal directions. 

Equation ^ with = is identical to the Dirac equa- 
tion for a free charged fermion in an s-wave background. 
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The Green's function for fermions in the condensed phase 
of a holographic s-wave geometry was studied by the au- 
thors of [26l - f28] for various values of the charge q and 
at zero temperature. The salient features of their anal- 
ysis are that the spectral function p has support inside 
an "infrared lightcone" which, in the probe limit, would 
correspond to the region cj > |fc|. In addition p may 
have non-vanishing support along codimension-one sur- 
faces outside this lightcone. The intersection of these 
codimension-one surfaces with the w = plane defines 
Fermi momenta kp. The location of the codimension- 
one surfaces coincide with the normal modes of Qi and 
C,2'- Suppose that at ky = 0, /c^, > 0, Ci has a normal mode 
which contributes to a pole of Gr,ii- Then, at ky = 0, 
k^ < 0, C2 will have a normal mode which contributes to a 
pole of Gr^22- As expected, the spectral function p is ro- 
tationally invariant. From ([t]) it should be clear that the 
locations of the normal modes of C, can be obtained from 
those of C with the replacement (w, kx) — ^ (— —kx). 

We would like to use the results of pSl - BS] as a starting 
point to study the distribution of normal modes of the 
Dirac equation ^ in the presence of a non trivial cou- 
pling 77. Unfortunately, the analysis of [29l [3^ implies 
that, in the probe approximation and with 7^ 0, one 
does not have a good zero temperature limit of the spin 
two condensate. Therefore, in what follows, we will re- 
strict ourselves to finite temperature configurations. At 
small but non-vanishing temperature, one expects the 
normal modes to broaden into quasi-normal modes whose 
width decreases with decreasing temperature. We have 
verified this expectation numerically. 

Consider first a configuration where 77 = and (^sa = 
and focus on a single normal mode of, say, Qi and the 
corresponding mode of C2. These modes cross at w = 
and kx = kp- Turning on 77 couples C,i and C2 and eigen- 
value repulsion replaces this degeneracy by a gap. This 
gap-generating feature of ([5| is similar to the mechanism 
described in |27j , where a gap was generated in an s-wave 
superconductor by adding a Majorana-like coupling to 
the fermion action. Next, consider a configuration where 
the condensate is rotated by 7r/4 [i.e. ipxy = 0). When 
77 > 0, Ci and Ci are coupled. Since their zero modes do 
not cross along the a; = plane a gap is not generated. 
For other angles one expects a continuous transition from 
the ungapped point (the node) to the gapped point (the 
anti-node). Generically, multiple FSs will be observed for 
a positive mass spin two condensate. When 77 these 
multiple FSs will appear as an intricate angular depen- 
dent band-like structure. In what follows we focus on a 
single FS. An analysis of the lifting of the degeneracy in 
the presence of multiple FSs is left for future work. 

Results. - We solved ^ numerically using a shoot- 
ing algorithm, first to generate the condensate and then 
to compute the Green's function. We set i^m^ = 7/4 
(corresponding to a d-wave order parameter of conformal 
dimension 7/2), 777 = and 77 = 0.02 qL. Numerically we 



found that once 77 becomes too large, there exists a crit- 
ical temperature below which the spectral function ex- 
hibits a gap even at 6* = 7r/4. In this sense, the value of 1] 
was tuned to be small. In presenting our results, we set 
ipA — and let 9 be the angle between k and the a:-axis. 

In fig. [I] we plot the spectral function p{uj, k) along 
the a; = plane. As expected at low temperatures, the 
gap-generating mechanism discussed above reduces the 
surface of quasi-normal modes intersecting the uj — 
plane to four nodal points. As depicted in fig. |4j close 
to the nodes and at small \oj/pq\ > 0, the quasi-normal 
modes have the structure of asymmetric Dirac cones. At 
the node, one can characterize the dispersion relation of 
the quasi-normal modes by the Fermi velocities v± = 
and V|| = Here is the momentum 

in the direction perpendicular (parallel) to the FS at the 
node, and uj{k) specifies the location of the quasi-normal 
mode. We find that v±/v\\ ~ 13. 

To compute the size of the gap along the FS as a func- 
tion of angle and temperature, we need first to identify 
the FS. Since we are working at non-zero temperature 
and since the spectral function is gapped, the FS is some- 
what ill-defined. We define the Fermi momentum kp{9) 
as the momentum at which the gap is minimized, or as 
the momentum at which the spectral function is maxi- 
mized if a gap does not exist. This definition is similar 
to the one in e.g. (TS]. We find that kp is almost inde- 
pendent of 9. A curious feature of our system is that kp 
decreases (by roughly 20%) as the temperature is lowered 
approaching kp ^ 0.11 g/i at low temperatures. 

The spectral function p{uj,kp{9), 9) is called the en- 
ergy distribution curve (EDC). In fig. [2]we plot the EDO 
for various angles, and for three different temperatures. 
We define the gap A{9,T) as half the energy between 
the two maxima of the EDC. By definition, A = when 
there is only one maximum. At temperatures above 
Tare — 0-56 ± 0.01 Tc thermal effects broaden the gap 
so that it is unobservable at all angles (fig. [2^). For 
temperatures T < Tgap ~ 0.28 ± 0.01 Tc (and with an an- 
gular resolution of 1/100) a gap is observed everywhere 
except for the node located at = 7r/4 (fig. [2]:). In this 
case, A{9, T) = Ao(T)| cos(2^)| to very good accuracy as 
expected from mean field theory [37) . At intermediate 
temperatures the gap opens up only outside an interval 
\0~l\<9o[T) (fig.[2]3). The interval where the spectral 
function remains ungapped is called a Fermi arc |14H16) . 
A plot of the dependence of the gap A on angle and tem- 
perature is shown in fig. [3j The dependance of the arc 
length on temperature is given in fig. [2] 

Discussion. - It is interesting to compare the holo- 
graphic model discussed here with (ARPES) photoemis- 
sion experiments on high Tc superconductors [3T1I33) . 
Such a comparison is problematic if one wants to get 
quantitative insight into the superconductive phase of the 
cuprates: Our holographic setup is merely a toy model 
for a true holographic d-wave superconductor with a well 
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FIG. 2: (Color online) EDCs, p{lo, kpiO), 6), for several angles 
and temperatures. In the left panel, figures (a), (b) and (c) 
correspond to temperatures T/Tc = 0.66, 0.43, 0.27 respec- 
tively. The bottom right plot shows the dependence of the 
Fermi arc length on temperature. The FS at T = 0.43 Tc is 
plotted on the top right. 
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FIG. 3: (Color online) Plots 
of the gap A(6l,r) at dif- 
ferent temperatures and an- 
gles. When the FS is com- 
pletely gapped, A can be fit 
by Ao(r)cos26l (dashed). 
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FIG. 4: (Color onhne) Max- 
ima of the spectral function p 
at T = O.lSTc along the indi- 
cated constant lu slices. The 
ratio of the Fermi velocities is 
wx/w|| — 13. The dashed line 
specifies the FS. 



defined boundary theory; the fermions described here are 
not, a priori, related to the fermions forming the conden- 
sate; there is no underlying lattice in our model. Nev- 
ertheless it is thought-provoking that, even in the ab- 
sence of several significant features of the cuprates, the 
holographic model manages to reproduce some of their 
notable properties. 

One key feature of cuprate superconductors is the ap- 
pearance of a d-wave gap in the superconductive phase 
. The gap A is well fit by the function | cos kx — cos ky \ 
evaluated on the FS, as expected from a square lattice, 
and is essentially indistinguishable from the mean field 
result A = Aocos(20). The four nodes associated with 
such a gap and the Dirac cones emanating from them 
have been observed experimentally These features 

are similar to the ones we observe in our holographic 
setup, depicted in fig. [3] and |4j A typical experimental 
value of the ratio between the quasi-particle Fermi veloc- 
ities is v±/v^^ ~ 20 [32]. Such a ratio can be obtained in 
the holographic setup by reducing ry. 



Fermi arcs are observed in the pseudogap phase of high 
Tc materials, at temperatures Tc < T < T* |14h16) . As 
described above, these are finite angular intervals along 
which the FS remains ungapped. The holographic Fermi 
arcs displayed in fig. [T] and |2] are distinct in that they 
appear in the superconductive phase, at temperatures 
Tgap < T < Tg.^c < Tc- In addition, the dependence of 
the arc length on the temperature (fig. [2] lower right) is 
certainly nonlinear whereas in [T5] it was observed that 
the arc length depends linearly on the temperature. 
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